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Abstract. Thomson's theorem states that static charge distributions in conductors 
only exist at the conducting surfaces in an cquipotcntial configuration, yielding a 
minimal electrostatic energy. In this work we present a proof for this theorem based on 
the variational principle. Furthermore, an analogue statement for magnetic systems 
is also introduced and proven: the stored magnetic field energy reaches the minimum 
value for superficial current distributions so that the magnetic vector potential points 
in the same direction as the surface current. This is the counterpart of Thomson's 
theorem for the magnetic field. The result agrees with the fact that currents in 
superconductors are confined near the surface and indicates that the distinction 
between superconductors and hypothetical perfect conductors is fictitious. 
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1. Introduction 



According to Thomson's theorem [TJ, given a certain number of conductors, each one 
with a given charge, the charges distribute themselves on the conductor surfaces in order 
to minimize the electrostatic energy. Thomson's result has been applied several times 
in the literature. In particular, in papers of this journal it was used 'to find the density 
of the induced surface charge in the classical case of a point charge in front of an infinite 
planar conductor' [10] and how 'the use of the minimum energy principle to explain 
static distributions of electric charge provides interesting physical insights' [3J. Other 
interesting applications of the Thomson's theorem are also found in [21 H] . Even though 
the proof of the theorem may be found with great detail in El [7] , a different approach 
not found in the literature, based on the variational principle, is presented in this paper. 

Moreover, the analogue situation for the magnetic field is also introduced: The 
magnetic field energy is minimal for superficial currents distributions so that the 
magnetic vector potential points in the same direction as the surface current. However, 
since energy conservation is assumed, this result cannot be applied to ohmic conductors 
in which a steady current can only exist if there is an electric field. Thus, it can only 
be used in classical perfect conductors and superconductors. Furthermore, it must be 
stressed that this theorem is not incompatible with Greiner's and Essen's result on the 
magnetic field energy maximization [HI [9] . Their work concerns currents in wires that 
can move in space while this study regards fixed currents. We shall prove this statement 
and Thomson's theorem by using a variational principle for the electromagnetic field 
energy functional. 



2. The proof of the theorem 

2.1. Thomson's theorem 

In equilibrium, the static energy functional for the electric field in a system composed 
of a conducting region surrounded by vacuum, may be written as: 

U = ^{ P 0-|(V0) 2 }d^ + J^dS 

= L W - ! (^) 2 1 dv - L ! (^) 2 dv + L ^ dS 

= J [p<P - | [V • (0V0) - 0V 2 0] } dV 

- f ~ { v • Uv<f>) - <pv 2 <p) dV+ f o-^dS (1) 

The volume of the conductor is represented by V in , the outside volume is V^, ut and S is 
the conductor boundary surface. Inside the conductor, we consider the volume charge 
distribution p and on the conductor's surface, the sufarce charge distribution a. After 
applying Green's Theorem, the energy functional becomes: 
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where V + and V~ are the gradient operators at the surface in the upper and lower 
limits, respectively. Since the total charge in the conductor is limited and constrained 
to the conductor volume and surface, one must include this restriction by introducing 
a Lagrange multiplier. Therefore, the infinitesimal variation of the energy reads: 

5<f)e V 2 (j)dV 



5U 



Vout 



{5<t> (p + e o V 2 0) + 5p (0 - A)} dV 

+ f [5(f) [a + e n ■ (v + - V - ^)] + 5a ((f) - A)} dS 
From the energy minimization, the Euler-Lagrange equations are: 



s 



A 



a = -e n ■ (v + - V~0) 
= A 



V out :V 2 







(3) 



(4) 



(5) 



(6) 



From these equations, in the minimum energy state, the potential is constant 
inside the conductor which means the electric charge distribution is superficial in an 
equipotential configuration proving Thomson's theorem. 

2.2. Minimal Magnetic Energy theorem 

The same procedure will be applied to the magnetic field. The energy functional for the 
magnetic field with no time dependencies is written as: 

= [ {]■ A- — (V x A) 2 \dV- I —(VxA) 2 dV+[k-AdS 
Jv in { 2// v / J J Vout 2/i v / J s 



V in I 2/iQ 



Ax (V x A)] +A - [V x (V x 1)]) \dV 



+ / k-AdS 
Js 



(7) 



Again, the volume of the conductor portion is represented by Vj n , the outside volume 
is V^t an d S is the conductor boundary surface. Inside the conductor, we consider 
the volume current distribution j and on the conductor's surface, the sufarce current 

— * 

distribution k. Using Green's Theorem, the energy functional becomes: 



{A- 




V x Vxi 



}dV 
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As in the electric field case, some constraints should be imposed. In particular, due 
to charge conservation, the electric current must obey the continuity equation both in 
volume and at the surface [TT] : 

V-J =0 (9) 

V s -k = (10) 

where Vs is the surface gradient operator. Notice that constraints are local, not global 
ones. In other words, the Lagrange multipliers are not constant but space functions. 
Therefore, the infinitesimal variation of the magnetic energy reads: 




+ 5j ■ [A - VA dV 



^5k- [A- V S AJ jdS(ll 

Hence, the Euler-Lagrange equations for this energy functional are: 
V x (V x A) = V x B = /i j 




x A - V~ x A 



V x B = 



(12) 

(13) 
(14) 



From eq. (12), the magnetic field must be zero inside the conductor, discarding 
any currents in volume (j = 0) in the minimum energy state. On the other hand, from 
eq. (13), one concludes that the potential vector A is tangent to the conducting surface. 



Since the magnetic field is zero in the inner volume, considering a cartesian coordinate 

system for a given surface element, if one of the axis is orthogonal to the surface element, 

the surface current becomes: 
1 „ 



k 



-n x 



v.. 





1 A 




) x V S A 


= — n x 


Vn x V S A 





(15) 



where Vs points in the same direction of A = V^A by definition and V5/ is orthogonal. 
From this result, it is clear that the surface current k will point in the same direction of 
A. Notice that no restriction was imposed on the vector potential. Therefore this result 
is independent of the chosen gauge. Even though the proof of the theorem does not 
include the presence of an external magnetic field, the result can be easily generalized. 
The presence of an external magnetic field demands another surface integral in eq ^ 
for currents at infinity. It does not affect the final result, however. 
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3. Applications 

3.1. Coaxial cable 

The Minimal Magnetic Energy Theorem can be illustrated in a simple system such as 
the coaxial cable. A coaxial cable can be modeled by an outer cylindrical conducting 
surface with radius b, carrying an electric current /, and a massive cylindrical conductor 
with radius a, carrying the same electric current in the opposite direction. Furthermore, 
for this case study, the inner current is divided into a superficial current I\ at r = a 
and a current I2 in the inner volume, so that I — Ii + I 2 . The magnetic field is readily 
obtained using Ampere's Law: 

w>h r 2 , 



B 



2<+> < r < a 



2nr a 

a<r<b (16) 
b < r 



Thus, the magnetic field energy on a length L of the cable is given by: 



U 




-L/"(^) 2Lnrdr + ^[ b (^-) 2 
2/i Jo \2iir a 2 J 2fi J a \2nrJ 



2Lirrdr 



= ^r + ^r log6/a (17) 

The magnetic energy reaches the minimum for J 2 = and therefore, for I = I\. In 
other words, the magnetic energy attains its minimum value when the current only flows 
at the surface. Furthermore, the current density was considered to be homogeneous 
in the azimuthal coordinate, which is intuitively the minimum energy configuration, 
leading to a tangent magnetic vector potential. A similar analysis for an electrostatic 
distribution can be found in 

3.2. Superconducting cylinder in a constant magnetic field 

The Minimal Magnetic Energy Theorem can also be useful to determine the 
configuration of the magnetic field and surface current distribution of a superconductor 
in equilibrium, for example, a cylindrical perfect conductor/superconductor with radius 
R immersed in a external constant magnetic field (see Figure 1). In order to achieve the 
magnetic vector potential solution, one must solve the following differential equations: 

VxB = VxfVxl)=0 (18) 



However, instead of starting solving immideatly, one should look into the 
symmetries of the system. For instance, if the external constant magnetic field points 
in the same direction of the y-axis and if the cylinder is considered to be infinite, there 
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Figure 1. Infinite superconducting cylinder in a constant magnetic field 



won't be any dependencies on the z coordinate (cylinder axis). Therefore, the magnetic 
field simplifies to: 



p d(j) dp 

Moreover, due to the symmetry of the system, the z-component of the magnetic 
field must be zero: 

^)-W =0 (20) 

These assumptions and constraints transform eq. (18) into: 

VxB- k( ld2Az 92 Az ldA A <) ( 2i) 
\p 2 d(f) 2 dp 2 p dp J 

which is the well known Laplace equation. Before writing down the general solution, let's 
analyze the boundary conditions. At p — > 00, the magnetic field must be equal to the 
external one: B cxt = B j = B ((j) cos (p + psin0). Furthermore, since the magnetic field 
is zero inside the perfect conductor, from eq. (19) one concludes the magnetic potential 
vector must be constant inside the perfect conductor/superconductor. Therefore, the 
solution for this case study is then: 

A z = const. + B ( p ] cos0 (22) 

VP / 

Even though this is only the solution for the z-component of A, the other components 
are not relevant to compute the magnetic field. However, the two remaining solutions 
can always be determined assuming the Coulomb gauge and using eq. (20) while the 
boundary condition comes from eq. ( 15 ) since the surface current must point in the same 
direction of the potential vector. Now everything is ready to calculate the magnetic field 
outside the cylinder: 

B = pB f 1 — ^\ sin 4> + 4>B f 1 + ^\ cos (23) 
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which implies that, 

B(p = R) = 2B cos <p } (24) 

The knowledge of the magnetic field on the cylinder's surface allows determination of 
the surface current: 

k = 2 — - cos k (25) 

The total current obtained through integration of the surface current is zero as expected, 
otherwise the energy would diverge. A more detailed analysis on this problem can be 
found in 



3. 3. Superconducting sphere in a constant magnetic field 

The same calculations can be performed for a superconducting sphere with radius R in 
a constant magnetic field pointing in the z-axis direction (see Figure 2). 
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Figure 2. Superconducting sphere in a constant magnetic field 



Similarly to the cylinder case, eq. (18) is simplified a lot by implementing the 
symmetries of the system. Since the external constant magnetic field points in the 
same direction of the z-axis, there won't be any dependencies in the <fi coordinate and 
the magnetic field along this coordinate must be zero. Therefore, the magnetic field 
simplifies to: 

' <} (^sin^f-^CVOfl (26) 



B 



r sin 8 36 



r dr 



Again, these assumptions and constraints transform eq. QlSp into: 

d 



V x B 




(Af sin I 



(27) 



Since the magnetic field is zero inside the sphere, from eq (26), the magnetic vector 
potential must have the following form: 

C 



r sin 9 



(2* 
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where C is a constant. In order to prevent the magnetic vector potential from diverging 
at r = and 6 = 0, the constant C must be zero. Furthermore, at r — > oo, the magnetic 
field must be equal to the external one: _B ex t = B$k = Bo(fcos6 — ^sin^). Therefore, 
the solution of eq (|27|) for this case is: 



A^r>R) = ^^r-^\wn9 (29) 
which leads to the following magnetic field outside the sphere, 

B = rB (l - ^ cos9 - 9B (l + ^) sinfl (30) 
From this result, the magnetic field at the sphere surface becomes: 

B = --Bowa.e§ (31) 



Using eq (13), the surface current density is therefore: 

£=- 35 °sin#0 (32) 
2 fa 

Unlike the infinite superconducting cylinder in a constant magnetic field, the 
superconducting sphere has a total non-zero electric current, / = 3^R. This current 
is essential in order to keep the magnetic field from entering inside the spherical 
superconductor. A similar approach to this problem can also be found in |13j . 



4. Conclusions 



In conclusion, the energy approach to the nature of magnetic fields has provided several 
interesting physical results such as [HI [151 EE EH HBj- in particular our theorem, 
analogously to Thomson's theorem in electrostatics, allows the prediction of surface 
current distributions in superconductors at equilibrium. Furthermore, it may change 
the distinction between superconductors and perfect conductors. 

In the pedagogical literature on superconductivity a distinction between 
superconductors and hypothetical perfect conductors is sometimes made. Both have 
zero resistivity below some critical temperature but when cooled below this temperature, 
while penetrated by constant external magnetic fields, the superconductor expels the 
field (the experimentally observed Meissner effect) but, in the corresponding thought 
experiment, the perfect conductor does not (see e.g. Poole et al. [19]). From the Minimal 
Magnetic Energy Theorem, one concludes that energy minimization and expelling 
external fields are the same thing. Therefore, the perfect conductors would in fact 
behave like superconductors in their approach to thermodynamic equilibrium. 
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